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ABSTRACT 

With the aim of studying active region fan loops using observations from the Hinode EUV Imaging 
Spectrometer (EIS) and Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA), 
we investigate a number of inconsistencies in modeling the absolute intensities of Fe VIII and Si VII 
lines, and address why spectroheliograms formed from these lines look very similar despite the fact 
that ionization equilibrium calculations suggest that they have significantly different formation tem- 
peratures: log(Te/i^) = 5.6 and 5.8, respectively. These issues are important to resolve because 
confidence has been undermined in their use for differential emission measure (DEM) analysis, and 
Fe VIII is the main contributor to the AIA 131 A channel at low temperatures. Furthermore, the 
strong Fe VIII 185.213 A and Si VII 275.368 A lines are the best EIS lines to use for velocity studies in 
the transition region, and for assigning the correct temperature to velocity measurements in the fans. 
We find that the Fe VIII 185.213 A line is particularly sensitive to the slope of the DEM, leading to 
disproportionate changes in its effective formation temperature. If the DEM has a steep gradient in 
the \og{Te/K) = 5.6 to 5.8 temperature range, or is strongly peaked, Fe VIII 185.213 A and Si VII 
275.368 A will be formed at the same temperature. We show that this effect explains the similarity of 
these images in the fans. Furthermore, we show that the most recent ionization balance compilations 
resolve the discrepancies in absolute intensities. With these difficulties overcome, we combine EIS and 
AIA data to determine the temperature structure of a number of fan loops and find that they have 
peak temperatures of 0.8-1.2MK. The EIS data indicate that the temperature distribution has a finite 
(but narrow) width < \og{(TT^/ K) = 5.5 which, in one detailed case, is found to broaden substan- 
tially towards the loop base. AIA and EIS yield similar results on the temperature, emission measure 
magnitude, and thermal distribution in the fans, though sometimes the AIA data suggest a relatively 
larger thermal width. The result is that both the Fe VIII 185.213 A and Si VII 275.368 A lines are 
formed at log (Te/K) ~ 5.9 in the fans, and the AIA 131 A response also shifts to this temperature. 
Subject headings: Sun: corona — Sun: UV radiation — Techniques: spectroscopic 



1. INTRODUCTION 

To understand how the solar corona is heated to high 
temperatures, it is important to explain the heating of 
closed field structures such as active region loops. There 
have been extensive studies of these structures, and re- 
ce nt progress and o utsta nding issues have been reviewed 
bv lKlimchukI (|2QQ6[ ) and lReald (|2Q1Q[ ). A key diagnostic 
of the heating of coronal loops is the differential emis- 
sion measure (DEM) distribution, because many coronal 
heating models make specific predictions as to its form 
and shape. For example, nanofiare reconnection mod- 
els (Parkerl 119831 . 11988 ) predict the prese nce of a weak 
high temperature component in the DEM (^Cargill|[l995.) . 
There have been sever al recent studies that have tried 
to detect this emission (jSchmelz et al.|[2QQ9l : iReale et"al] 
[M)9l : [TbstiTenillMc^ 

The DEM gradient, or proportion of hot and cool 
material, also sets constraints on i mpulsive, steady or 
quasi-steady loop heating models () Warren et al.l I2Q1QI : 
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iTripathi et al.|[2"QlQ[ ). and the nanofiare model also pre- 
dicts a spread in the temperature distribution within a 
loop due to the incoherent heating and cooling of unre- 
solved threads. There has been considerable debate as to 
whether loops have a multi-thermal temperature distri- 
bution because of conflicting measurements by different 
instruments (Lenz et"al] Il999l : lAschwanden et al.l 119991 : 
Schmelz et al. 200l|). The discuss ion hinges on issues 
such as background subtraction (jPel Zanna fc MasonI 
l2QQ3l ). methodo logv (lAschwandeiJ 120021) . or the spa- 
tial resolution (jAschwanden et al.l I2008D, or tempera- 
ture resolution of the instruments used ("Martens e t al.l 
l2002f) . Recent observations by the Hinode (Kosugi e t al.l 
12007) EUV Imaging Spectrometer (Culhane et al. M2007l . 
EIS) suggest that loops formed near IMK have a nar- 
row temperature d istribution, but are not isothermal 
(| Warren et al.l[2008h . 

Doppler velocity measurements are another important 
diagnostic of the heating process, and again there have 
been numerous recent studies of flows in active region 
loops using EIS data (Doschek et al. 2007; Hara et al] 
[20081 : iDel Zannal [2008: Brooks fc Warren 2009) . Some 
of the signatures of coronal heating models are expected 
to be quite subtle, however. The nanoflare reconnection 
model predicts weak downflows at 'warm' temperatures 
(jPatsourakos fc Klimchud l2006l ). and short-lived faint 
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upflows at high temperatures (iPatsourakos fc KhmchukI 
I2OO9). The signatures of nanofiare heating in spectral 
line profiles therefore depend sensitively on the temper- 
ature of the flows, and this also implies that accurate 
measurements will allow inference of the properties of 
the energy release. It is important therefore to assign 
the flow temperatures as accurately as possible. 

A class of loops that have not yet been stud- 
ied in sufficient detail with the latest instrumenta- 
tion are the fan structures that appear as partially 
observed long loops at the edges of active regions. 
They are s een mostly in the 0.4-1. 3M K tempera- 
ture ran^e (|Schriiver et a l.l 1 19991 : iDel Zanna fc MasonI 
I2QQ3I : lUgarte-Urra et a l. 2009) and have densities greater 
than l oR {NeJcm-^) 9 (Del Zannal 12003: 'Young et al." 
l2007bD. Th ev also appear to sho w red- shifted downfiows 
(jWinebarger et al.ll2Q02l : iMarsch e Tall [20041 ). We exam- 
ine a sample of these structures in this paper. 

For this and all the above observational studies, spec- 
tral emission from ions of Fe is of major importance. The 
high elemental abundance of Fe leads to many strong 
emission lines in spectrometer and imager pass-bands, 
yet interpretation of the observations is often difficult 
because of uncertainties in th e atomic data used in pre- 
diction of emission (L anzafame et al.l 120021 : 1 Young et al.l 
[2009). In principle, the emission from spectral lines of 
Fe can provide stringent constraints on temperatures and 
densities in the corona, however, it is of paramount im- 
portance that the diagnostic capabilities of Fe lines be 
assessed critically. 

Initial results from EIS have indicated a number of 
problem s in interpreting the Fe emission. I Young et al.l 
(|2007b[ ) noted that spectroheliograms of Fe VIII and 
Si VII look very similar in specific active region features 
such as the fan loops, despite the fact that the tempera- 
tures of the peak fractional abundance in ionization equi- 
librium are significantly differen t; \o^{TeJ K) = 5.6 and 
5.8, respectively, according to iMazzotta et al.l ([1998). 
Since Fe is a considerably more complex atom than Si 
and is therefore likely to be less well understood, this 
observation has led to the suggestion that the ionization 
balance calculations for Fe V III need to be revis ed up- 
ward to higher temperatures (|Young et al.ll200"7b[ ). Fur- 
thermore, this could have an impact on surrou nding ions 
such as Fe VII and Fe IX. Recently, Young ~fc Landil 
(|2009f ) found evidence that this is indeed the case for 
Fe VII. 

In addition to this temperature problem, a number 
of previous differential emission measure (DEM) studies 
have found difficulties reproducing the absolute inten- 
sities of Fe VIII and Si VII lines simultaneously. See, 
for example, the q uiet Sun off-limb DEM analysis of 
IWarren & Brooks (20091). the on disk quiet Sun study 
of Brooks et al.. (.20091). or the analysis of a cool active 
region feature bv iLandi fc YoungI ()2009f ). 

Such issues are important to resolve as the Fe VIII 
185.213 A and Si VII 275.368 A lines would provide the 
best EIS lower temperature constraints on the DEM if 
we had good confidence in the atomic data. Also, since 
they are strong, they are the best lines to use for velocity 
measurements in the transition region, and are present 
in the majority of EIS observations. As discussed, it is 
therefore of critical importance that we have confidence 



in the temperatures we assign to the measured velocities. 

Changes to the formation temperatures or atomic data 
for these ions could also have an impact on interpreting 
the results from imagers. For example, changes to Fe VIII 
could affect the response functions for the 131 A channel 
of the Atmospheric Imaging Assembly (AIA) on the Solar 
Dynamics Observatory (SDO). 

Motivated by our interest in studying the temperature 
and velocity structure of the fan loops, this situation has 
led us to take a closer look at the formation of the Fe VIII 
and Si VII spectral lines. Recently, there have been a 
number of revisions to the ionization balance calculations 
(jBryans et al.ll2009l : iDere et ani2009f ) and we investigate 
whether they could help resolve these inconsistencies. We 
also examine the details of spectral line formation in this 
temperature range in the quiet Sun and in the fan loops. 
In particular, we examine whether convolving the contri- 
bution functions with the temperature distribution of the 
feature could explain the similarity of Fe VIII and Si VII 
images. It is known that the temperature of peak contri- 
bution to the line intensity can be shifted from the theo- 
retical peak temperature of the emissivit y if the shape of 
the DEM is taken into consideration (Brosius e t al.ll l996l : 
iFeldman et al .111 9991 : iDel Zanna et ani2003f ). In doine; so, 
we finally show that the ionization equilibrium calcula- 
tions for Fe may not be the source of this problem. By 
determining more realistic effective formation tempera- 
tures, we present a possible explanation for the obser- 
vations. We also show that the most recent ionization 
balance calculations can resolve the discrepancies in the 
magnitudes of the intensities found in previous studies. 

With confidence in the atomic data for these ions re- 
stored, we perform an emission measure (EM) analysis of 
a number of fan loops using EIS and AIA and compare 
the results from the two instruments. 

2. SIMILARITY OF FE VIII AND SI VII IMAGES 
Figure [T] shows example images of AR 10978 in (top 
left to bottom right) Mg V 276.579 A, Mg VI 268.986 A, 
Mg VII 280.737 A, Si VII 275.368 A, Fe VIII 194.663 A, 
Fe IX 197.862 A, and Fe X 184.536 A. These images 
are formed at \og{Te/K) = 5.4, 5.6 5.8, 5.8, 5.6 
5.8, a nd 6.0, respectively, according to IMazzotta et al.l 
(|1998r ). The data were obtained on 2007, December 10, 
at 00:19:27UT. The EIS instrument has 4 slit options 
(l'^ 2'^ 40'^ and 266^0- The V' slit was used for these 
observations and stepped over a FOV of 460'^ by 384'^ 
The exposure time was 40s at each position. The data 
were processed using standard EIS data reduction rou- 
tines available in SolarSoft (eis_prep). 

It is clear that the images do not follow the expected 
temperature trend. The Fe VIII and Si VII images look 
very similar, especially in the fan structures to solar east 
of the AR and the cool loops to solar west. This sug- 
gests a similar temperature of formation. In addition, 
the Mg VI and Fe IX images look to be formed at lower 
and higher temperatures, respectively. Note, for exam- 
ple, the third column of the figure. The vertical extent of 
the cool features to solar west and the fan structures to 
solar east appears to increase from top to bottom (Mg VI 
to Si VII to Fe IX) suggesting an increase in tempe ra- 
tures. The problem identified bv I Young et al.l (|2007b[ ) is 
clearly seen in these examples. 
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Fig. 1. — Hinode EIS rasters of AR 10978 showing the similarity of Fe VIII and Si VII images. The Mg VI and Fe IX images look to 
be formed at lower and higher temperatures, respectively. For example, note the increasing vertical extent of the fan structures and cool 
features to east and west of the active region as you look at the images from top to bottom of the third column (Mg VI, Si VII, Fe IX). The 
images are scaled linearly. 



3. ATOMIC DATA 

In this work we use several sources of atomic data 
as follows. All the contribution functions are calcu- 
lated using the CHIANTI v6.0.1 datab ase (iDere et al. 
I1997L 120091 ). The coronal abundances of lFeldman et al. 
(Il992f) were used, together with the ionization fractions 
of ' Mazzotta et alJ (|1998l ) and lDere etld] (j2009f ). 

We use four specific lines to explore how their for- 
mation temperatures change in the quiet Sun: Mg VI 
268.986 A, Fe VIII 185.213 A, Si VII 275.368 A, and 
Fe IX 188.485 A. For these lines we u sed the effective 
collisi on strengths and radiative data of Landi & Bhatia 
(I2007D. [Griffin et a l. (2000), Bhatia & Landi (2003), and 
iStorev et al.l ()2002f ). respectively. A large number of ad- 
ditional spectral lines were used to compute DEMs and 
the data sources for these secondary data are too nu- 
merous to list here. We refer the interested reader to 
the CHIANTI database where the references are given in 
detail. 



4. SPECTRAL LINE FORMATION 
TEMPERATURES 

The intensity of an optically thin emission line arising 
from an atomic transition from level i to level j can be 
expressed as 

hj = j G,j{Te,Ne)cl){Te)dTe (1) 



where the differential emission measure (DEM) is defined 
as, 

(/)(Te) = N^dz/dTe (2) 

with dz the differential distance along the line of sight 
and A^e the electron density. The contribution function 
is defined as 

G,j{Te,Ne) = A{Z)F{Te,Ne)e,j{Te,Ne)j^^ (3) 

with A{Z) the elemental abundance, uh the hydrogen 
number density, eij(Te^Ne) the line emissivity normal- 
ized to the ground state population, and F{Te^Ne) the 
ionization fraction as a function of temperature and den- 
sity. The accuracy of the computation of F for Fe and 
Si i s important for this p aper. In Figure [2] we contrast 
the lMazzotta et ahl ()1998f ) ionization equilbrium calcula- 
tio ns for Fe VIII an d Si VII with the updated results 
of iDere et aTl (j2009f ). Note that the peak formation 
temperature does not change, so the temperature prob- 
lem remains. The maximum fractional abundance, how- 
ever, increases by 10% for Si VII and 55% for Fe VIII. 
These changes could be important when modeling the 
magnitudes of the absolute intensities. Note that the 
Bryans et al. (2009) calculations give the same forma- 
tion t emperatures for Si VII and Fe VIII as De re et al.l 
pool . The maximum fractional abundances are also 
within 10%. 

The function G is strongly peaked at the characteris- 
tic temperature of the line, but the temperature of peak 
contribution to the total intensity can be shifted from 
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Fig. 2. — Comparison of t he Fe VIII and S i VII i onization equi- 
librium cur ves calculated bv lMazzotta et all ([l993) with those of 
IDere et al.l (^009). 



the theoretical one if there is relatively more emitting 
material at higher or lower temperatures. Another key 
issue for this paper is that if the contribution function is 
broad, there may be significant contributions to the total 
intensity from multiple temperatures. In Figure [21 note 
that the Fe VIII cur ve is signific antly broader than that 
of Si VII in the new iDere et all (|2009) calculation, and 
encompasses much of the temperature range of that ion. 
A further refinement to the analysis then, is to calculate 
the effective formation temperature (Tg//) of the line by 
convolving the contribution function with the DEM dis- 
tribution for the target region or structure of interest, 
and then taking the temperature at which that function 
has its maximum value. 



^(Te//,7Ve) = max(^(Te,7Ve)) 
Vg) is obtained by differentiat 

giTe,Ne) = GijiTe,NemTe) 



(4) 

where g{Te,Ne) is obtained by differentiating Equation 

n viz. 

(5) 



4.1. Formation Temperatures in the Quiet Sun 

Before turning to the AR fan loops, we first investi- 
gate the formation temperatures of lines in the quiet 
Sun. Since we are also interested in how the more re- 
cent ionization balance calculations affect the derived 
DEM, and whether they could help in resolving known 
discrepancies, we derived a new QS DEM using the lat- 
est updated atomic data from v6.0.1 of the CHIANTI 
dat abase and a dopting the ionization fractions calculated 
by iDere et ahl (|2009). Throughout this paper we use 
the pre- la unch laboratory photometric calibration of EIS 
(jLang et a l. 2006). 

Perhaps the best way to test the accuracy of the atomic 
data is to use off-limb QS observations where there is less 
contamination along the line of sight from multiple tem- 
peratures. Therefore, we selected lines for analysis that 
were identified as yielding consistent emission measure 



distributions above the q uiet solar limb. Most of them 
are taken from Table 4 of IWarren Brooksl (j2009[ ). The 
critical issue here, however, is the consistency of lines in 
the temperature range where Fe VIII, Fe IX, and Si VII 
are formed, and it would be best if we had independent 
constraints on the DEM that do not involve these lines. 
Therefore, we used deep exposure observations of an on- 
disk QS dataset. Deep exposure data are helpful for 
bringing out the weaker cool lines below log(Te/i^) = 
5.8, that would be harder to detect off-limb, and there- 
fore for constraining the DEM at the lower temperature 
end. This also allows us to directly compare the results 
to that of lBrooks et al.l ()2009t) where various inconsisten- 
cies in observed and predicted intensities were found in 
on-disk data using the older ionization balance calcula- 
tions. 

The data used were taken on 2009 June 14 starting 
at 09:59:42UT. Context images of the quiet region are 
shown in Figure [3] including both Si VII 275.368 A and 
Fe VIII 185.213 A, which again look very similar. The 
EIS study used the 2" slit with 4 min exposures and 
scanned a FOV of 60''x 512^' in about 2hrs 10 mins. 
About 70% of the wavelength range of the CCDs were 
read-out and telemetered to the ground. Since this 
does not cover th e full C CD, some of the lines used in 
IWarren fc Brooksl (|2009f ) were not available. In addi- 
tion, we included several important lines of interest for 
this work, i.e., Mg VI 268.986 A, Fe VIII 185.213 A, and 
Si VII 275.368 A, and, as discussed, additional cooler 
temperature lines to provide constraints on the DEM 
(O IV 279.631 A, O IV 279.933 A, O V 248.456 A, Fe VIII 
186.601 A, Fe VIII 194.663 A, Si VII 272.641 A, Mg VII 
278.402 A, and Mg VII 280.737 A). Ah the lines used are 
listed in Table [T] which also shows which lines were se- 
lected for the analysis of the fan loops in Sections [5] and 

m 

The data were processed using standard procedures 
available in Solar Soft. In addition, however, the orbital 
variation of the line centroids and spectral line tilt were 
estimated by single Gaussian fits to the Fe XII 195.119 A 
line and removed from the data prior to further analysis. 
The offsets in X- and Y- between the short and long 
wavelength CCDs were corrected for so that the same 
area was averaged. Thes e offsets were estim ated using 
the software developed by "Young et al." (2009). The data 
were then averaged over the common area of the two 
CCDs and fitted using single or multiple Gaussians as 
appropriate. 

4.2. Differential Emission Measure 
The method for determination of the DEM distri- 



bution is the sa m e as described in Warren (2005^, 
iBrooks fc WarrenI (|2006[ ). and Brooks et al. (2009^ . 
Briefiy, we represent the DEM with a series of spline 
knots, the magnitudes of which are found by min- 
imization of the differences between the observed and 
DEM predicted intensities. This minimization is per- 
forme d by the SolarSoft routine MPFIT (M arkwardtl 
l2009l ). The number and initial values of the spline knots 
can be pre-set and they can be interactively manipulated 
to control the smoothness of the emission measure. This 
technique has been found to be useful for representing 
rapid changes in the DEM. 
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Fig. 3.— Top row: Context observations of the quiet Sun region. From left to right, images in Fe VIII 185.213 A, Si VII 275.368 A, Fe X 
184.536 A, Fe XI 188.216 A, Fe XII 195.119 A, Fe XIII 202.044 A, Fe XIV 274.203 A, and Fe XV 284.160 A. The images are scaled linearly 
between their minimum and maximum values, except for Fe XV which is scaled logarithmically because it is very weak. Images longward 
of 250 A have been shifted upward to account for the offsets between the EIS CCDs. The common area of the CCDs was averaged to 
compute the DEM. The Y-dimension has been rebinned for presentation. Lower row: Left panel, comparison between the quiet Sun DEM 
derived using the latest atomic data in this paper and that of Brooks et al. ( 2009). Right panel, normalized g(Te, Ne) function in the quiet 
Sun for Mg VI 268.986 A, Fe VIII 185.213 A, Si VII 275.368 A, and Fe IX 188.485 A. The lines are formed at log (Te/K) = 5.68, 5.75, 5.85, 
and 6.05, respectively. 



We measured the electron density using the following 
three line ratios: Fe XIII 196.525 A/202.044 A, Fe XIII 
202.044A/203.826A, and Fe XII 186.880 A/195.119 A. 
There are a number of issues with blending and the 
accuracy of the atomic data within certain ranges for 
these diagno s tic ra tios that are discussed in detail by 
I Young et al.l (|2009f ). None of these issues are expected 
to affect our measurements significantly, however. The 
ratios give densities for the quiet Sun of log {Ne/cm~^) 



= 8.7-9.0, indicating a pressure of log {Pe/cm~^K~^) 
= 14.9-15.2. The upper ends of these ranges are from 
Fe X II, which tends to give higher values Young et al.l 
|2009[ ). Since we cannot rule out that the possibility that 
this is a real density change with temperature, we choose 
to assume that the electron pressure is constant. One of 
these assumptions is necessary to cast the intensity inte- 
gral in the form of Equation [1] and calculate the contribu- 
tion functions. We adopt a value of log {Pe/cm~^K~^) = 
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15.0 in this work and compute the contribution functions 
for ah the emission hues of Table [T] usin g CHI ANTL 

The quiet Sun DEM of iBrooks et~all (j2QQ9f ) was used 
to give a first guess for the temperature knots and spline 
fit to the QS data. The fit was then adjusted slightly 
to improve the agreement between observed and pre- 
dicted intensities. The results are given in Table [T] and 
shown in comparison with our previous result in Fig- 
ure [31 80% of the lines are reproduced to within 30%. 
Fe XIV 270.519 A and Fe XV 284.160 A are significantly 
overestimated but these lines are very weak in the spec- 
trum and the DEM is falling rapidly in this temperature 
range. Fe XIII 200.021 A, Fe XIV 264.787 A, and Fe XIII 
197.434 A are underestimated by the DEM but these re- 
sults are consistent with our previous work, and the latter 
is blended with Fe VIII 197.362 A 



(jBrown et al.M2008 ). 
Mg VII 278.402 A is underestimated by about a factor of 
two, but this discrepancy can be resolved if the contri- 
bution from th e Si VII 278.445 A bl end is accounted for 
as described in I Young et al.l ()2007bl ) . 

All of the Fe VIII, Fe IX, and Si VII lines are con- 
sistent with each other and the predicted intensities 
are within 30% of the observed ones. These results 
are interesting because the DEM curve derived using 
the updated atomic data and ionization fractions from 
CHIANTI v6.0.1 resolves the inconsistencies found us- 
ing older versions in several previous studies of Fe VIII 
and Si VII lines in the quiet S un off limb, on disk, and 
in cool active re^on features (fWarren fc Brooks 2009 : 
IBrooks et aH 120091 : iLandi fc Youndl2009f ). 

4.3. Discussion 
Using Equations |4] and [5] we determined Tg/ / in the 
quiet Sun for four key lines (Mg VI 268.986 A, Fe VIII 
185.213 A, Si VII 275.368 A, and Fe IX 188.485 A). The 
gradient of the DEM slope in the quiet Sun increases in 
the temperature range of formation of these lines (Fig- 
ure [3]) so that there is relatively more material above 
the temperature of the maximum ion abundance for all 
of them. This results in all of these lines being formed 
at temperatures above their maximum ion abundance. 
As noted, however, the Fe VIII 185.213 A line is par- 
ticulary sensitive to this because of the larger width of 
the F(Te, Ne) curve. It is therefore disproportionately af- 
fected, and is actually formed much closer in temperature 
to Si VII 275.368 A than the ionization balance would 
suggest. An illustrative example is shown in Figure [3l 
Note that the g{Te,Ne) curve for Fe VIII 185.213 A is 
considerably broader than for the other lines, suggesting 
that there will be contributions from multiple tempera- 
tures. Note also that because the contribution functions 
for individual lines differ within an ion they may also 
have different effective temperatures. 

5. FORMATION TEMPERATURES IN ACTIVE 
REGION COOL FAN STRUCTURES 

In the previous section we discussed how the slope of 
the DEM could lead to the Fe VIII 185.213 A and Si VII 
275.368 A lines being formed at similar temperatures in 
the quiet Sun. This DEM analysis was made, however, 
on averaged spectra over a relatively large area, and the 
differences may not be dramatic because the structure of 



the atmosphere does not change rapidly over the small 
temperature interval (log {T^/ K) = 5.75-5.85) where the 
lines are formed. It is clear, however, that specific fea- 
tures such as the cool fan loops look very similar in Si VII 
and Fe VIII images. It remains to be seen whether this 
explanation could hold for these specific features of in- 
terest. 

To address this question we investigate the tem- 
perature distribution of the extended fan loops in 
AR 10978 (shown in Figure [1]). This active re- 
gion has p revio usly been studied by several authors 
(iDoschek et al.H2008l: [Brooks et al.l 120081: IWarren et al.l 
20081: lUgarte-Urra et all 120091 : iBrvans et al.l 120101 : 
Brooks fc Warrenll2011f ). 



5.1. Context images and intensity ratio map 

Context images of the region in both Fe VIII 185.213 A, 
and Si VII 275.368 A are shown in Figure HI As with the 
analysis of the QS, we corrected for the CCD offsets, 
orbital variation and spectral line tilt prior to further 
analysis. In this case we use the artificial neural net- 
work model of Kamio et al.l ((2010f ) to correct the orbital 
variation of the line centroids and grating tilt. The rest 
of the data processing was accomplished using standard 
procedures (eis.prep). The fan region is shown as a large 
box in Figure IH 

Also shown in Figure 2] is a map of the Fe VIII 
185.213 A/Si VII 275.368 A intensity ratio. This quanti- 
tatively describes the similarity of the two images. Note 
that the ratio appears to vary in the active region core. 
Fe VIII 185.213 A is blended with N i XVI 185.251 A 
(|Young et al.ll2QQ7at [Brown et al.l l2008V which is formed 
at 2.8MK, and it likely contributes to this emission and 
influences the ratio. In contrast, one can see that the fan 
region is almost uniformly dark, indicating little varia- 
tion in the ratio across this area. A histogram of the 
ratio values for the pixels within the boxed region is also 
shown in the Figure. A Gaussian fit to this distribution 
indicates that the mean value is 2.3 and the standard de- 
viation is 0.34, i.e., the variation in the ratio across the 
fans is less than 15%. 

5.2. EM analysis 

We computed the emission measure in several small 
boxes extending along one of the fan loops (shown as re- 
gions E1-E6 in Figure [4]). The intensities were averaged 
in these areas. The background emission was estimated 
for each box by averaging the intensities in the adja- 
cent boxes also shown in the Figure. Single and multiple 
Gaussian fits were made as appropriate. 

For this analysis we used most of the same lines as 
we used for the QS. These are indicated in Table [H A 
few lines were not available in the current dataset and 
could not be included. These lines can also be iden- 
tified from the Table. Furthermore, we included O VI 
183.937 A, O VI 184.117 A, Fe XVI 262.984 A, and Fe XVI 
265.003 A to provide additional low and high tempera- 
t ure constraints . 

I Warren et al.' ("2008") analyzed a sample of 'warm' EUV 
loops in this region by fitting a Gaussian distribution in 
temperature for the model emission measure of the form. 



EMo 



exp 



(6) 
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Si VII 275.352A 



Fe VIII 185.2 




Fig. 4. — Example showing the similarity of Fe VIII 185.213 A and Si VII 275.368 A images. A ratio map and histogram are shown 
indicating the low variation of the ratio across the fans. The histogram is made from the region indicated by the large box. An enlargment 
of the boxed region is also shown to highlight the locations used for the analysis along the fan loop (E1-E6). The adjacent boxes show the 
areas used for background subtraction. 



For comparison with this work we use the same fitting 
technique here. The fit is made by minimization of the 
differences between the measured and computed intensi- 
ties. The density is ahowed to be a free parameter in the 
fitting and we obtain values of log {Ne/cm~^) = 9.2-9.7, 
decreasing from base to apex. These results are compa- 
rable to the d ecre ase along another example reported by 
lYoung et al.l pOOTbh . 

The computed temperature distributions for positions 
E1-E6 along the loop together with their peak tem- 
peratures, emission measure, and Gaussian width are 
shown in the top two rows of Figure [5l There are a 
few problematic lines that are not well reproduced by 
the model. They vary by position, but generally include 
O IV 279.933 A, O V 248.456 A, and O VI 183.937 A, ah 
of which are underestimated to varying degrees. Fur- 
thermore, the background subtracted high temperature 
lines of Fe XIV and Fe XVI are often underestimated, 
including the normally robust Fe XIII 202.044 A line. 
O V 248.456 A is blended with an Al VIII line that is 



a weak contributor in most conditions, hence the agree- 
ment found for O V in the QS. The blend, however, could 
contribute more when the emisison measure is peaked 
mostly at coronal temperatures. The Mg VII 278.402 A 
and 280.737 A lines also sometimes show unusual behav- 
ior. In most cases they are reproduced well, but mis- 
matches of 50-80% are sometimes seen. In these cases, 
the discrepancies for Mg VII 278.402 A cannot be re- 
solved by accounting for the blend with Si VII 278.445 A. 
Since it is mostly the low and high temperature lines that 
show discrepancies it could be that the temperature dis- 
tribution deviates from a Gaussian shape in the wings. 
In general most of the lines are reproduced well, however. 
In total, about 80% of the lines are reproduced to within 
40% with most better than 30%. 

5.3. Discussion 

Compared to the quiet Sun, the fans show a significant 
enhancement of emission measure (factor of 3-7). The 
distribution peaks at log (Te/K) = 6.0 near the base (El) 
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Fig. 5. — Emission measure results along the loop from the base (top left) to the top (bottom right) and the effect on the intensities 
of Fe VIII 185.213 A and Si VII 275.368 A. Top two rows: temperature distribution with the peak temperature, emission measure, and 
Gaussian width indicated. Bottom two rows: Intensity of the two spectral lines as a function of temperature. The vertical dotted lines 
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and increases along the loop to a maximum of log {Te/K) 
= 6.1 (E3 and E4) before decreasing again. The dis- 
persion in the Gaussian distribution appears to decrease 
along the loop length. Previous studies of cool loops in 
active regions have suggested that they are isothermal 
fPel _Zanna & Mason 2003). Our results indicate that at 
its narrowest, the width of the emission measure distri- 
bution is still log (cTTe/i^) ~ 5.4. In co niparison with 
the ^w arm' loops in this region studied by I Warren et al.l 
(|2QQ8l ). this fan loop has a broader temperature distribu- 
tion along most of its length. Only towards the tip (E6) 
does it approach the values for the 'warm' loops. Note 
that the gradient of the slope of the emission measure 
distribution in the temperature range log {Te/K) = 5.6- 
5.8 is greater than in the quiet Sun, even at the base of 
the fan (El). 

Figure [5] also shows the normalized giTe^Ne) function 
for the Fe VIII 185.213 A, and Si VII 275.368 A lines 
that results from convolving their contribution functions 
with the fan emission measure at each position. We see 
that even when the temperature distribution is relatively 
broad (El), the effective formation temperatures of the 
lines are separated by only 0.1 in the log. It also indicates 
that the separation in temperatures reduces as we go 
along the fan loop until they are both formed at the same 
temperature of \og{Teff / K) = 5.85. Clearly when the 
EM distribution is narrow in this temperature range the 
two lines will be formed close together in temperature. 
This analysis also shows that this is true even if the EM 
distribution is relatively broad, provided it has a positive 
gradient. The Fe VIII 185.213 A line seems to be very 
sensitive to the slope of the EM distribution in general. 

6. EMISSION MEASURE ANALYSIS WITH EIS 
AND AIA 

To examine whether our explanation could hold in gen- 
eral for all fan loops we investigated a number of other 
cases. The loops we selected were observed in active 
region 11093 on 2010, August 13. In thi s analysis, w e 
also include data obtained by SDO/AIA (|Golubl [20061 ). 
As mentioned in the introduction, it is also important 
to understand the formation of spectral lines in detail 
if we are to correctly interpret observations by EUV 
imagers. The broad wavelength pass-bands of such in- 
struments include multiple lines and their contributions 
may change in quiet and active conditions. We illus- 
trate this here by examining the temperature response 
of the AIA 131 A filter. The pass-band of this filter is 
expected to contain contributions from Fe VIII 130.941 A 
and 131.240 A, Fe XX 132.850 A, Fe XXI 128.755 A, and 
Fe XXIII 132.906 A with the Fe VIII lines the dominant 
contributor outside of flares. As we have seen, the for- 
mation of Fe VIII lines is non-trivial. 

Figure [6] shows full Sun images from AIA taken at 
lOUT. AR 11093 was the target for Hinode and the EIS 
FOV is overlaid as a box (120'' by 512''). The EIS ob- 
servation started at 09:26:40 and lasted about 1 hour. 
During this time the 1" slit was moved in steps of 2" 
across the region with 60s exposures at each position. 
The EIS data were reduced and processed in the same 
way as discussed previously. 

The AIA data were obtained pre-processed from the 
level 1.5 cutout service. These data have had flat-fielding 



and dark current corrections applied. Bad pixels and 
cosmic ray hits have also been repaired by an iterative 
method. We converted the data to DN/pixel/s. The im- 
ages have also been rotated to solar North and shifted 
to match the SDO/HMI (Helioseismic and Magnetic Im- 
ager) FOV center so that they are all coaligned. Nev- 
ertheless, we detected some small offsets between dif- 
ferent filters and corrected this misalignment by cross- 
correlating the full disk images. 

Coalignment of the EIS and AIA data was achieved 
by extracting the EIS FOV from the AIA 131 A fuh Sun 
image and cross-correlating it with the Fe VIII 185.213 A 
EIS raster. Due to the orbital variation of the posi- 
tion of the EIS slit, the spatial sampling is not uniform 
across the raster. In contrast, the AIA image is taken 
nearly instantaneously (2.9s exposure) so there is negli- 
gible movement. Therefore, the effective plate-scale mag- 
nification between the images is not uniform in the so- 
lar X-direction. For these observations this discrepancy 
amounted to 1" across half the raster and was corrected 
manually. 

Figure [7] shows images from EIS and AIA of the EIS 
FOV. For this analysis we selected small boxes on each of 
3 loops and averaged the intensities in these areas. The 
selected positions are shown in the figure as EAl, EA2, 
and EA3. The background emission was estimated for 
each box by averaging the emission in the adjacent boxes 
(also shown in the figure). Single and multiple Gaussian 
fits to the EIS data were made as appropriate. The spec- 
tral line-list for these observations was smaller than that 
of the December 2007 region, but it still covered a suf- 
ficiently broad temperature range from log (Tg/i^) = 5.5 
to 6.4. The selected lines are indicated in Table [H In 
addition, we included O VI 184.117A, Mg V 276.579 A, 
and Mg VI 270.394 A to provide further lower tempera- 
ture constraints. 

The emission measure analysis is the same as that used 
in the previous section for consistency. Since the back- 
ground subtraction sometimes reduces the intensities of 
the density diagnostic lines to zero, we have assumed a 
density of \og{Ne/cm~^) = 9.5 in computing the emis- 
sion measure for both EIS and AIA. This value is the 
average of that obtained along the loop analyzed in §5.21 

We extracted intensities from the coaligned 131 A, 
171 A, 193 A, 211 A, 335 A, and 94 A filters. These filters 
are expected to be dominated by fines of Fe VIII, Fe IX, 
Fe XII, Fe XIV, Fe XVI, and Fe XVIII. The 94 A filter 
also contains contributions from Fe VIII and Fe X which 
can dominate under certain conditions (jO'Dwver et al] 
2010). We use only the intensities that correlate well 
with the 131 A intensity for generating the emission mea- 
sure distribution. In practise, this generally excludes 
the 335 A and 94 A bands where there is little emis- 
sion for the fan loops we study. A photometric cali- 
bration uncertainty of 25% is used as the intensity er- 
ror (Boerner et al., in preparation). The response func- 
tions were calculated for a grid of electron densities 
and temperatures spanning log (A^e/cm~^) = 6-12 and 
log {Te/K) = 4-8. This was done by calculating isother- 
mal spectra as a function of wavelength at each temper- 
ature and density, ^(A, Te, A^e), using CHIANT I v6.0. 1 
and the coronal abundances of iFeldman et all (|1992l ). 
The spectra were then convolved with the effective ar- 
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AIA 131 A 13-Aug-2010 10:00:09 




AIA 193 A 13-Aug-2010 10:00:07 



Fig. 6. Full Sun images from SDO in the AIA 133 A, 171 A, 193 

El 

eas for each filter, £^i(A), and integrated in wavelength, 
J Ei{X)S{X^Te^ Ne)dX^ to obtain the response functions, 
Ri(Te^Ne). The effective areas were obtained from So- 
larSoft. 

The EM distributions computed from the EIS data for 
EAl, EA2, and EAS are shown in the top row of Figure 
[HI The results computed from the AIA data are shown 
in the second row. EM loci curves are also shown to 
indicate the sensitivity of the EIS lines and AIA filters. 
Table [2] shows the peak temperature (logTo), emission 
measure (log£^Mo), and Gaussian width (logao) derived 
from EIS and AIA for each loop. The uncertainties in 
these parameters are also given. 

The two instruments agree that the peak temperatures 
of the loops are log (Tq/ K) = 5.9-6.0. The emission mea- 
sure magnitudes are \og{EM/cm-^) = 26-27 and the 
EIS and AIA results are consistent to 25-55%. The EIS 
data indicate that one of the loops (EAl) is isothermal 
and that the other two have a finite (but narrow) width. 




AIA 171 A 13-Aug-2010 10:00:12 




, and 211 A filters. The overlaid box shows the FOV used in Figure 

consistent with the findings for the loop in §5.2[ The 
AIA data yield similar results. For the two loops that 
are not isothermal, AIA finds a slightly larger thermal 
width (0.14 in log). This may be due to the coarser tem- 
perature resolution compared to EIS. The uncertainties 
in all parameters from AIA are also larger. Note that an 
error of zero is returned in the isothermal case because 
the temperature becomes fixed. 

Again in Figure [8] we show the intensities of Fe VIII 
185.213 A and Si VII 275.368 A as a function of tem- 
perature after convolution with the EM distribution of 
each loop. In all three cases the two lines are effectively 
formed at the same temperature. 

We also investigated the effective response of the 
AIA 131 A filter after convolution with the EIS emis- 
sion measure distribution. The main contributing lines 
in this pass band at log (Te/K) < 6.4 are Fe VIII 
130.941 A, and Fe VIII 131.240 A. These produce the 
peak around log (Tg/i^) = 5.7 in the usual response. 
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Fig. 7. AIA and EIS images of fan loops. From the left: AIA 133 A, EIS Fe VIII 185.213 A, AIA 171 A, EIS Fe X 184.536 A, AIA 193 A, 
EIS Fe XII 195.119 A, AIA 211 A filters, and EIS Fe XIV 264.787 A. The images are scaled logarithmically. 



^131 (^e) (dashed line in Figure [8|). The convolved func- 
tion Ri3i {Te)^{Te)Te is plotted in Figure [HI as the solid 
line. Note that the peak of the response shifts to higher 
temperatures after being convolved with the EM distri- 
bution; log(Te/i^) = 5.85-5.95. This is consistent with 
the analysis in the previous sections and demonstrates 
that these effects can also be important for imager ob- 
servations. 

7. SUMMARY 

Motivated by the desire to study the temperature 
structure of active region fan loops we have attempted to 
resolve inconsistencies found in previous work using EIS 
data. In particular, we have shown that the similarity 
in EIS Fe VIII 185.213 A and Si VII 275.368 A images, 
that is not expected from the respective temperatures of 
peak abundance in ionization equilibrium, can be under- 
stood when a more accurate calculation of the effective 
formation temperature in the solar corona is performed. 
This is done by convolving the contribution functions 
with the DEM of the target of interest. If the DEM has 
a steep gradient in the log (Te/K) = 5.6-5.8 range, or 
is sharply peaked, the two lines will be formed close in 
temperature. In this work, we compared the effect of this 
technique on the formation temperatures of these lines 
in the quiet Sun. The initial separation of log (Te/K) 
= 5.6-5.8 is reduced to 5.75-5.85, and it is clear that 
Fe VIII 185.213 A has a substantial contribution from 
emitting material at log {Te/K) = 5.8. 

To examine whether this explanation could work for 
the fan loops, we derived the EM distribution along one 
example in AR 10978. The temperature distribution 
peaks near 1-1. 2MK and narrows along the loop. The 
peak contribution to the line intensity is log {Te/K) = 5.9 
for both the Fe VIII 185.213 A and Si VII 275.368 A lines. 
To investigate whether this effect is generally applicable 



to other AR fans, we examined a number of other loops. 
We found that in all cases the two lines are formed at 
the same temperature (log Te/K ~ 5.9). This suggests 
therefore, that the expected difference in images of the 
fans formed from these lines is a result of an overestima- 
tion of the separation in formation temperatures by the 
approximate method of assuming the lines are formed 
at the temperature of the peak fractional abundance in 
ionization equilibrium. 

Note that other lines may be affected in similar ways. 
For example, in Figure [T] the Mg VII and Fe IX images 
look different despite the fact that they have similar ion- 
ization equilibrium temperatures. We have verified that 
they are in fact formed at different temperatures in the 
fan loop of Section [5l 

To demonstrate the importance of understanding the 
formation of the EUV spectrum for broad pass-band im- 
agers we studied the effect of convolving the AIA 131 A 
response function with our EIS fan loop EM distribu- 
tions. We showed that, as a result, the peak of the dom- 
inant lower temperature part of the effective response 
function shifts up to log {Te/K) ~ 5.9. 

It is important to emphasize that we have not shown 
that this explanation holds in general for all active ar- 
eas or structures on the Sun. If the DEM slope is shal- 
lower (or flat) in the log {Te/K) = 5.6-5.8 range then 
the Fe VIII 185.213 A and Si VII 275.368 A lines should 
still be formed at a wider separation in temperatures 
and examples of significantly different images should be 
found. The apparent lack of such observations for any 
solar feature, however, provides a stringent constraint 
on the gradient of the DEM slope in this temperature 
range all over the Sun. This is consistent with other 
independent studies that show strong similarities in the 
shape of th e DEM distribution in d ifferent areas of the 
quiet Sun (jLanzafame et al.l I2QQ5I : iBrooks et all I2QQ9I : 
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Fig. 8. — Emission measure results for the three loops EAl, EA2, and EA3 from EIS and AIA and the effect on the intensities of Fe VIII 
185.213 A, Si VII 275.368 A, and the AIA 131 A response function. Top row: EIS EM distributions. Second row: AIA EM distributions. 
EM loci curves are overlaid in color. Third row: Intensity of the two EIS spectral lines as a function of temperature. The vertical dotted 
lines indicate the formation temperatures. The separation in formation temperatures is shown in the legend. Bottom row: The AIA 131 A 
effective response as a function of temperature. The original response is shown by the dashed line. 
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iFeldman etaH [20091 ). 

The DEM-gradient resolves the inconsistencies in for- 
mation temperatures between the Fe VIII and Si VII 
hues, but a related issue is that discrepancies have also 
been found in the magnitude of the intensities of these 
lines in previous DEM studies using EIS data. We 
showed here that these additional issues are resolved 
when t he most re c ent io nization balance compilation 
data of iDere et al.l ()2009l ) are used for the atomic cal- 
culations. From this analysis, therefore, no substantive 
relative error in the Fe ionization balance is indicated for 
the specific ions of Fe VIII and Fe IX. 

This work has demonstrated that the strong Fe VIII 
185.213 A and Si VII 275.368 A lines can be used with 
confidence for DEM studies and velocity work in the 
transition region. Therefore, we examined the temper- 
ature structure of a small sample of fan loops. We 
found that they have peak temperatures in the range 0.8- 
1.2MK. One loop was found to be isothermal, but more 
often the temperature distribution has a narr ow width. 
This result is similar to that found by War ren et al.l 
(|2008) for 'warm' active region loops. In one detailed 
case, the EM distribution is found to broaden consider- 
ably towards the base. This could have implications for 
the location of the heating. 

We also found that the peak temperatures and emis- 
sion measures derived from AIA data are in agreement 



with those derived from EIS. There is also agreement 
on whether the loops are isothermal or not. The AIA 
analysis indicates a slightly larger thermal width than 
EIS when the loops are not isothermal. This is possibly 
because the EIS data contain observations from consec- 
utive ionization stages of Fe whereas the AIA data only 
sample every second ionization stage. 
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TABLE 1 

Quiet Sun intensities and line selection. 
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Intensities are reported for lines used in the QS DEM analysis. The ticks 
indicate the lines that were used for the QS DEM (Eq^), and the EIS (Ejan) 
and EIS-AIA {EAf an) analyses of the fans in sections [S] and (6] 

Units are erg cm~^ s"-*^ sr"-*^. 



TABLE 2 

Comparison of DEM results from EIS and AIA. 







EIS 






AIA 




Loop 


log To 


log EMq 


log CTQ 


log To 


log EMq 


log cro 


EA 1 
EA 2 
EA 3 


5.93±0.01 
5.90±0.04 
6.05±0.01 


27.10±0.02 
25.88±0.09 
26.44±0.02 


4.50±0.00 
5.41±0.18 
5.46±0.04 


5.88±0.05 
5.97±0.13 
6.03±0.08 


27.18±0.07 
26.10±0.10 
26.33±0.10 


4.50±0.00 
5.55±0.20 
5.60±0.27 



